STABILITY OF SOLUTIONS OF QUASILINEAR 
PARABOLIC EQUATIONS 

GIUSEPPE MARIA COCLITE AND HELGE HOLDEN 

Abstract. We bound the difference between solutions u and v of ut = aAu + 
div x f + h and vt = bAv + div x g + k with initial data if and ij>, respectively, 
by \\u(t,-)-v{t,-)\\ LP(E) <A B (t)Wp-<l>\\l P £ (R n) +B(t)(||o-6||«, +||V«-/- 
Vi ■ g\\ao + \\f u - g u \\oo + \\h — fc||oo) Pp |E|' )p . Here all functions a, f, and h 
are smooth and bounded, and may depend on u, x £ 1", and t. The functions 
a and h may in addition depend on Vti. Identical assumptions hold for the 
functions that determine the solutions v. Furthermore, E C R n is assumed to 
be a bounded set, and p p and r) p are fractions that depend on n and p. The 
diffusion coefficients a and b are assumed to be strictly positive and the initial 
data are smooth. 



1. Introduction 

We show that one can bound the difference between solutions u and v of 

u t = a(t, x, u, Vu)Aw + diva; (/(£, x, u)) + h(t, x, u, Vti), x eR n , <t <T, 
u(0,x) = <p(x), x£R n , (1.1) 

and 

v t = b(t, x, v, \7v)Av + div x (g(t, x, v)) + k(t, x, v, Vv), x £ R", < t < T, 
v(0,x)=ip(x), xel", (1.2) 

respectively. The assumptions are that the diffusion coefficients a and b are bounded 
from below by a strictly positive constant. All functions a, /, h, etc, as well as the 
initial data ip, etc, are assumed to be smooth and bounded. We are interested in 
estimating the local L p -norm of u(t, •) — v(t, ■) over any bounded subset E C K™ in 
terms of norm differences of the initial data as well as a and b, etc. 

In the hyperbolic case, that is, a = b = 0, the classical result of Kuznetsov [*H2) 
and Lucier ( see a ls° [HI Ch. 2]) reads 

||tt(*, •) - v(t, -)|| L1(R) < \\ V - i>\\ LHM) + t min{T.V.(^), T.V.(^)} ||/ - 9|| Lip 

in the one-dimensional case (n — 1) where / = f(u), g — g(u) and h = k = 0. 
Here T.V.((j>) denotes the total variation of the function <f> and ||/|| Li p denotes the 
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Lipschitz semi-norm. Recently, Bianchini and Colombo 2 showed flux stability in 
the case of hyperbolic systems on the line. Indeed, they established the estimate 

\\u(t, ■) - v(t, .)|| L1(R) <Ct\\Df- Dg\\ CO{Q) 

for solutions u and v of ut + f(u) x = 0, Vt + g{u) x — 0, respectively with u\t=o — 
v\t=o- The usual assumptions on the flux functions and the initial conditions apply, 
see PJ. 

The dependence in a of the solution u of the equation 

u t — Ao(u) = 

is treated in pQ, assuming only that a is nondecreasing, and thereby allowing de- 
generate diffusion. However, no explicit stability estimate is provided. Otto |T5) 
studied the equation 

B(u) t - div K (a(Vu, B(u))) + h{B{u)) = 

with a continuous and monotone nondecreasing B. Under certain assumptions he 
proved that 

||B(ui(i)) - B{u 2 (t))\W < exp(Lt) ||B(«i(0)) - ^MO))^ . 

By extending Kruzkov's famous doubling of variables method, Bouchut and Perthame 
[3| showed that 

,o „,0|| 



\ui(t, ■ 



< M 



■CT.V.(u?)>/tLip(a) 



u 2(^, OIIl^r") - iri ~ "2ii L i 

when Uj satisfies ut + div x (/) = Aa(uj) with initial data Uj, j = 1,2. Here a is 
assumed to be Lipschitz and nondecreasing. 

Closer to the approach of this paper, Cockburn and Gripenberg established 
the estimate 

| il(Rn) <T.V.(v# \\fi-ti\ 



AVtn 



||ui(t,.)-U2(V)ll-i 

for solutions Uj, j — 1,2 of 

u jt t = div x (fj) + A(a,j(uj)), Uj\ t =Q = p. 

Allowing for explicit spatial dependence in the flux function, Evje, Karlsen, and 
Risebro showed stability for solutions of 

Uj, t + div x (kj(x)fj(u)) = AAj(u), Uj\ t =o = u°, 

in the sense that 



||tti(t,-)-ua(*,-)ll 



LH 



< w 



+ *C'(ll*i-Moo,6 + ll/i-/alloo,up) 



ViC 



A 1 y & 



Where II ' lloo.fo and 



oo Lip ^ s ^ ne sum °f t ne sup-norm and the BV-norm and the 
sum of sup-norm and the Lipschitz norm, respectively. Here Aj is allowed to be 
degenerate. Karlsen and Ohlberger established L 1 contractivity of solutions of 

u t + div a (V(t, x)f(u)) = V • (K(t, x)VA(u))) + q(t, x, u). 

Recently, Chen and Karlsen [S] established the estimate 



\\u 1 (t,-)-u 2 (t,-)\\ 



Elixir 



+tc\\f[-m 



1/2 
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for solutions of Uj t t + div x fj(iij) — V • (Aj(uj)Vuj) with initial data Uj\t—o = Uj. 

We consider here the strictly parabolic case where the diffusion constant is not 
allowed to decrease to zero. However, we allow full explicit spatial and temporal 
dependence in all parameters. In addition, we let the diffusion and source depend 
explicitly on the gradient of the unknown u. All parameters, including the initial 
data are assumed to be smooth. Existence of regular bounded solutions is secured 
by classical results, see The question is to obtain explicit stability estimates. 
Our main result reads as follows. Let u and v denote solutions of ljl.ll) and l|1.2[) . 
respectively. Then 

\\u(t, ■) - v(t, -)\\ LP{E) < A B {t)\\<p - ^ll^ (R „) 

+ B{t)(\\a - b\\ L ^ {Ko) + IIV* ■ / - V* • <?|U» (TC ) 

)Pp 
\E\*> , 



where 

A E (t) 



B{t) := C 



ifl<P<2, . = f^ + ^> ifl<P<2, 

if2<p<oo, p ' [i, if2<p<oo, 

|^|(2- P )/2p+l/2„ + | S |l/p^ if 1 <p<2, 

(1 + t (P-a)/P)( l^l 1 /^ + lE] 1 ^), if 2 < p < oo, 

t, if 1 < p < 2, 

(t + t 2 /P), if2<p<oo, 



for any bounded connected set E C K™ with Lipschitz boundary. Here TZq = 
[0,T] x Ex [-K u Ki] x [-K 2 ,K 2 ] and TL = [0,T] x E x [-K^K^. 
As a particular example we note that for solutions u and v of 

u t = a(t, x, u, Vu)Au, vt = b(t, x, v, Vu) Av 

with initial conditions u\t=o — and v\t=o — ip, we find 

||u(t, •) - v(t, -)\\ LHE) < cm 1 ' 2 - + iei 1 ' 2 ) y - tf|| £oo(RB) 

+ Ct|^||a-6||£ (Wo) . 

Our proof is based on a homotopy argument, inspired by 0j. Introducing 

u e , t = (9a + (1 - 9)b) Au e + div x (Of + (I - 9)g) + 9h + (1 - 9)k, 
u e \ t =n = 6<p+(l-6)<ip, 

we see that uq — u and u\ — v. Thus iig interpolates between u (for = 0) and v 
(for = 1). The key estimate establishes that 

||tt(i, ■) - v(t, -)Wlp(e) = distiP( E )(u(i, •)>«(*> ■)) 
< length LP(B) (we(t, •)) 



d9, 

LP(E) 



and we establish ^-independent estimates for \\dug/d9\\. 
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2. Fundamental assumptions 



Fix T > 0. Let u = u(t, x) and v — v(t, x) be the bounded solution of the 
quasilinear initial value problem (see |13j) 

u t = a(t, x, u, Vu)Au + div-r (/(*, x, it)) + h(t, x, u, Vit), iel°, < t < T, 
u(Q,x) = <p(x), iff, (2.1) 



and 



v t = b(t, x, v, Wv)Av + div x (g(t, x, v)) + k(t, x, v, Vu), x S R", < t < T, 



(2.2) 



v(0,x) = ip(x), 

respectively. Here 

/=(/i,..,ttF(si,.,S„):RxI"x 

and 



div x (f(t,x,u)) = ^2 (/jlM.^Vu))^ 

3=1 

= £(fe+f».<) 

3=1 

= V x -f + f u - Vlt. 

Observe that V s • / is a scalar. The divergence operator div^ always acts on the 
spatial variables only. By V 9 a (similarly for b, h, and k) we denote the gradient of 
a with respect to the final n variables (where Vm usually sits). Our fundamental 
assumptions are 

(Hi) the viscous coefficients a and b are of class C 3 ([0,T] x R" x R x R n ) such 
that 



< a» < a(-, ■, •, •) < a* < oo, ||a||c3([o,T]xR™xRxR") < h, 



< K < &(•> ■>•>■) <b* <CO, ||b|| C 3([0,T]xK"XKXM 

for some positive constants a„, a*, 6*, 6*, fci; 
(H2) the convective terms / and g are of class C 3 ([0,T] 

x R" x R x 



(2.3) 



i n x M) and the 
such that for all 



source terms h and k are of class C 3 ([0, T] 
i, i, I € {1, ... , n} and any $ e • • . , /„, gi,...,g n ,h, k} the following 
quantities 

<9 3 $ 



<9$ 




<9 2 $ 








dxidxj 


L 












9m 




du 2 




3 



dxidxjdxi 



(2.4) 



<9 2 $ 



dxidu 



|Vg/l||z,oc, ||V g fc|| L » 



are all bounded by a positive constant fc 2 ; 
(W3) the initial data (/? and ^ are of class C 2 {R n ) such that 

IM|c 2 (R"), IIV'llc^R") < fc 3, 

for a positive constant &3. 



(2.5) 
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Lemma 2.1 (L°°-bounds on u and v). Fix T > 0. By ^2] there exist positive 
constants K\, K2, K3 such that 



MIl=°([0,T]xR"), |MU=°([0,T]xK") < K±, 



du 



dx,. 



L°°([0,T]xR~) 



dv 



d 2 u 



dxidxj 



L°°([0,T]xR") 



dxi 
d 2 v 



L°°([0,T]xR") 



< K 1 



(2.6) 



dxidxj 



L°°([0,T]xR") 



<K 3 , 



for all i, j G {1, ...,n} where K\, K2, K3 depend only on T, n, a*, a*, 6*, b* , fci, 
&2 , and /c3 . 

3. The homotopy argument 

Our approach is based on the following homotopy argument. Let < 9 < 1. 
The function ug interpolates between the functions u and v. More precisely, denote 
by ug the solution of the quasilinear initial value problem 

ug tt = (Qa(t,x,u g ,Vu e ) + (1 - 6)b(t,x,ug, Vu e ))Au e 

+ div x (0f(t, x, ug) + (1 - 0)g(t, x, u e )) 

+ 6h(t, x, ug,Vug) + (1 - 9)k(t, x, u e ,Vug), x e K", < t < T, 
ug{Q,x) = 6cp(x) + (1 - 6)ip(x), (3.1) 
Clearly 

MO = V, 1*1 = U. 

Indeed 

9^ug(t,-) 
is a curve joining v(t, •) and u(t, •), and 

INV) ~ v (t, -)Wlp(e) = dist LP ( E) (u{t, -),v(t, •)) < lengthy •)), (3.2) 

for each < t < T, E C 1™ measurable set and 1 < p < 00. 

Lemma 3.1 (L°°-bounds on Ug). By |13l Theorem V 8.1], there exist positive 
constants K\, K%, K3 depending only on T, n, a*, a*, 6*, b* , fex, fc 2 and fc 3 suc/i 
i/iai 



M e|U°°([o,T]xR") < Ki, 
<K 5 , 

L°°([0,T]xR") 



(3.3) 



d 2 Ug 



dxidxj 



< K 



6- 



L°°([0,T]xR") 

/or eac/i < 6 < 1 and i, j 6 {1, n}. 

Lemma 3.2 (Smoothness of ue). Assume {'Hi), (H2), and (H3). T/ie curiae 

6e [0,1] ^ug(t,-) e C7 2 (M") 
is 0/ c/ass C 1 . In particular, we infer 



length LP(is)l 
/or each < t < T and E C K.™ measurable set. 



Lp(E) 



dff. 



(3.4) 
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Proof. Consider the map 

F: V — > C°°(]0,r[xK n ) nC 2 ([0,T] x R n ), 

F(6,uj)(t,x) :— ^(t,x) - (6a(t, x,uj(t, x), Vu(t, x)) 

+ (1 - 9)b(t,x,uj(t,x),Vuj(t,x))^jAoj(t,x) 

- div x (6f(t, x, uj(t, x)) + (1 - 0)g(t, x, uj(t, x))) 

- (6h{t, x, u(t, x),Vu>(t, x)) + (1 - 9)k(t, x, uj(t, x), Vw(t, x))) , 

where 

V := {(0,w) e [0,1] x C°°(]0,T[xM n ) nC 2 ([0,T] x R n ) | w(0,-) = + (1-6>)V>}- 
From the definition of ug, 

F(9,u e ) = 0, < 6» < 1. (3.5) 

Observe that F is of class C 1 and 
OP 

— (9,uj) = (6(t,s,w,Vw) — o(t,a;,w,Va)))Au 

+ divj; 2;, w) — /(f, x, cj)) + fc(t, w, Vw) — x, w, Vu). 

To compute 



f)F F)F 

— { e^)[{e',z)] = —{e,u + ez) 



s=0 



we find 

F(6, uj + ez) = ^ + e^- - (9a(t, x,uj + ez, Vuj + eVz) 
at at 

+ (1 - 6)b{t, x,uj + ez, Vw + eVz)) (Aw + eAz) 

- div x (0/(i, x, uj + ez) + (1 - %(*, x, uj + ez)) 

- (6h(t, x,uj + ez, Vw + eVz) + (1 - 9)k(t, x,uj + ez, Vw + eVz)) , 
OF Q z 

— {9,uj + ez) = — - (6a(t, x, uj + ez, Vw + eVz) 

+ (1 - 9)b(t, x,uj + ez, Vw + eV z)) Az 
, 3a 



- (6-^(t, x, uj + ez, Vuj + eVz) 



db 

+ (1 - 6) — (t,x,uj + ez, Vuj + eVz))z(Auj + eAz) 
ouj 

— (6V q a(t, x, uj + ez, Vuj + eVz) 

+ (1 - 9)V q b(t, x,uj + ez, Vuj + eVz) ) • Vz (Auj + eAz) 



div,. (j6^(t,x,uj + ez) + (1 - 6)^(t. r - _-.)).; 



(d—(t, x,uj + ez, Vuj + eVz) 



dk \ 
+ (1 - 9) — (t,x,uj + ez,Vuj + eVz))z 

OUJ J 
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, x, u + ez, Vui + eVz) 
+ (1 - 9)V q k(t, x,uj + ez, Vuo + eVz)) ■ Vz. 



Thus 

— (0, w) [(^', z)] = - (6a{t, x, w, Vw) + (1 - 0)6(t, w, Vw)) Az 

>^(t,x,uj,Vuj) + (l-9) — (t,x,u),Vuj)) zAuj 
aui ouj J 

- (9V q a(t, x, uj, Vlu) + (1 - 9)V q b(t, x, uj, Voj)) ■ Vz Auj 

- dW x ((e^t, x, w) + (1 - 9)^-(t, x, uj)) z) 

0—(t,x,w,Vw) + (1 - 9) — (t,x,w, Vuj))z 

x, u>, Vu>) + (1 — 9)Vqk(t, x, uj, Vui) J • Vz, 

(6,uj),(d',z)€V. 

Observe that (6', z ) E T> satisfies the equation 

BF 

— (0,u,)[(0',z)]=C 
if and only if z is solution of the linear initial value problem 
Zt = (Oa(t, x, uj, Vuj) + (1 — 9)b(t, x, uj, Vw))Az 

+ (9a u (t,x,u,Vcj) + (1 - 9)b lAj (t,x,uj,Vuj))Aujz 

x, uj, Vuj) + (1 — 9)\/ q b(t, x, uj, Vw)) ■ Vz Aw 

+ div K ( (9 fa, {t, x,uj) + (1- 9)g w {t, x, uj)) z) 
+ (0h„ (t, x,uj,Vuj) + {l-9)k u} (t, x, uj, Vw))z 
+ (9V q h(t, x, uj, Vuj) + (1 - 9)V q k(t, x, uj, Vw)) • Vz + ((t, x), 
x G K", < t < T, 
z(fl,i) = «VW + (l-^(4 arGM". 



Since this problem is well-posed (see Theorem IV 5.1]), ^(9,u>) is invertible 



dF i 

By the implicit function theorem, the curve 9 \ — ► ug is of class C 1 and clearly (|3.4() 
holds. This concludes the proof. □ 



Differentiating equation l|3.1|l with respect to 9, we have 

= (9a(t, x, ug, Vug) + (1 - 9)b(t, x, u e ,Vu e ))A 

( 'da , \ . dug 

+ \9—[t,X,Ug,Vug) + (1 - 9) — (t,X,Ug,Vug) I All 8 -7^- 

{d9, 

a(t, X, Ug,Vug) — b(t, X, Ug, Vug)) Aug (3-6) 



+ (9V q a(t, x, ug,Vu e ) + (1 - 9)V q b{t, x, u e , Vug)) • V ^— ) Aug 
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+ ( —(t,x,ue) - -^{t,x,u e )j Vttfl + V x ■ f(t,x,u$) - V x • g(t,x,ug) 
^( t> x |Ufl) + (l-^( tia :^)).v(^ 
9W X ■ |£(t, x, u e ) + (1 - 0)V X • u,)) ^ 

+ (^(*, «») + (i - <? )^(*, ■ v we ^ 

+ h(t, x, u e ,\7u e ) ~ k(t, x, ug,Vue) 

( n dh . . . . <9fc , \ dug 

+ I 9—{t,X,U9, Vue) + (1 - 0) — (t,x,u e , Vu e ) I 

+ (0V 9 fc(t, au, u e , Vu e ) + (1 - 0)V g fc(t, i, u e , Vue)) ■ V f ^ ) 



Denoting 

Zg(t,X 



a(0,t,x) := 9a +{1- 6)b, 

(3(9, t,x) := (6V q a+(l-6)V q b)Au e + 6^ + (l-6)^L 
+ (9V q h+(l-9)V q k), 

1 (9,t,x) := fe^ + (l-6)^\ Au e 
a(6, t, x) := (a - b)Au e + (jL - ^ . Vu 

+ f-Vx-g + h-k, 

for each < 6 < 1, t > 0, a; € R n , there results 

^ = a(0, i, x) Az e + /3(0, t, a;) • Vz + j(6, t, x)z e + a{9, t, x), 

0< < 1, < < < T, xe R". (3.7) 

Moreover, observe that 

z g (0,x) = ip(x)-ip(x), < < 1, x £ R n . (3.8) 

Lemma 3.3 (L°°-bounds on a, (3, 7). From i/ie definition of a, l|2.3f) and (|3.3|l . 

we have 

0< a* <a( v , •)<<**, l|Va|| L oo < + # 5 +nir 6 ) (3.9) 

w/iere 

a* := min{a*, 6*}, a* := max{a*,6*}. 
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Moreover, from the definition of (3 and (12.4(1 , we infer 

\\0\\ L oo= sup ll&IUoo < K 7 , (3.10) 

j = \,...,n 

where 

K 7 := nkiK§ + 2k 2 . 
Finally, from the definition 0/7, (|2.3[1 . (I2.4|l and (|3.3[1 . we find 

hh-<K&, (3.11) 

where 

K 8 := nkiK & + (n + 1 + nK 5 )k 2 . 

Lemma 3.4 (L°°-bounds on zg). Assume (Hi), (Ti. 2 ), and (H3). There exists 
a positive constant C\ depending only on T , n, a*, a*, 0*, 6* , k%, k 2 , and k^ such 
that 

\\ze(t, -)IU-o(ii») <Cit+\\<p- V>IU- ( Rn), (3.12) 
/or eac/i < t < T and < < 1. 

Proof. To simplify the notation we let w denote the solution of (|3.7|l . that is, 

w t = aAw + j3 ■ Viu + 7W + a, w\t=o — w o- (3.13) 
Linearity implies that 

w — Wi + W2 

where w\ and w 2 solve 

101,4 = aAwi + [3 ■ Vioi +7101, ioi| t= o = io , 
u> 2 ,t = aAw 2 + /3 • Vw 2 + ^yw 2 + a, w 2 \ t =o = 0, 

respectively. We infer from p. 389] that 

Wl (t,x) = ( G{t,Q,x,i)w Q {Odi, 



w 2 (t,x) 



f f G(t,T,x,0o-(T,0d^dr, 

JO JR™ 

where G is the Green's function. For t E [0, T] for some fixed T positive we find 

\w 2 (t,x)\ ^ctiML,. 

Introduce z = w\ — wq which satisfies the equation for w 2 with a — aAwo + 
13 diva; wq + 7u>o . Thus 

|t«i(t,a:)| < \z(t,x)\ + \w (x)\ < \\w \\ 00 + Ct\\aAv} +P-'Vv3o+'yw \\ 00 . 

a 

Observe that in the previous lemma, the smoothness of the initial condition 
enters in a crucial way. With less regularity we get the familiar O^ 1 / 2 ) behavior 
near t = (see, e.g., [3 Sec. 4.4]). 
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4. Stability of quasilinear parabolic equations 

We begin with the following lemma. 

Lemma 4.1 (Poincare-type inequality). There exists a positive constant Aq, 
depending only on n, such that 



I \f\ 2 dx<A \B\ 2/n [ \Vf\ 2 dx + A \B\ 1/n [ \f\ 2 dx, 

JB JB JdB 



(4.1) 



for each f G C 2 (R n ) and B C M" bounded connected set with Lipschitz boundary. 
In the case n = 1 we mean 

[ \f\ 2 dx = \f(x )\ 2 , 

JdB 

for some xq G B. 

The proof of this lemma is more or less classical (see ^3 Theorem A. 9] and 
Lemma A. 2]) and the dependence of the coefficients on the measure of the domain 
is consequence of a standard rescaling argument. 

Now we prove the key estimate in the L 2 -norm for the map zg. 

Lemma 4.2 (Case p = 2: Energy estimate). Assume (Hi), (H2), and (H3). 
Then there exists a positive constant C'2 depending only on T, n, a*, a*, &*, &*, k\, hi, 
and &3 such that 

\\ze{t, OII^B) < C 2 (\E\ 1/2n + \E\ 1/2 )\\ V ~ 

+ C 2 1 \E\ 1/2n (||a - 6|Uoo (TCo) + ||V X • / - V* ■ sIU-cr) (4.2) 

\ 1/2 

+ Wfu - 9u\\l<*>cr.) + \\ h - k\\L°°(n ) j 

for each 0<t<T,0<8<\ and E C K™ bounded connected set with Lipschitz 
boundary. Here K = [Q,T] x E x [-Ifi, #1] x [-K 2 ,K 2 ] a«d 71 = [0,T] x £7 x 
:-K x .K{\. 

Proof. Let B C R™ be a ball and < t < T. Then by we find 
^ y ^(M)^ = y zgzg^dx (4.3) 

= / Azg efe + / zg(3-\7zgdx+ / 72^ rfx + / azgdx. 

JB JB JB JB 

Observe that, by (|3~TT|) . 

/ JZg dx < K s / zg dx, (4.4) 

JB 

and, by ftlOjl . 

/ zgp ■ \7zg dx < — { \(3\ 2 z 2 g dx + ^ / |Vz e | 2 dx (4.5) 
Jb "* Jb 4 Jb 

<M f z 2 dx + ^ f \Vzg\ 2 dx. 
a* Jb 4 Jb 

By Lemma EPI and <f!T5|) . 

/ azgdx< / |er| |z#| da; = / v|^[(vM \zg\) dx 

JB JB JB 
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<~J s W\dx+±J\a\z 2 e dx (4-6) 



<\\B\ \\a\\ L ~ {n) + C 2 \B\t 2 \\a\\ L ~ {u) 
+ \\<P~ ^lli«.( R «) \B\ \\<t\\l<*>(K) 

Ah 



<^\B\ h\\ L ~ { n)+Ci\B\1?\\a\\ L <~ m , 



where 

Moreover, by the divergence theorem we have 



Kq := 1 1 °'- 2 



/azgAzgdx = / azg(Vzg ■ v) dx — / (Vq • Vzg)zg dx — / a\Wzg\ 2 dx 
_? Job Jb Jb 

< I az g (V 'ze ■ v) dx + -— / |Va| 2 z^ete (4.7) 
Jas 2a* J B 

+ / | Vze| 2 c?x — a* / \\7zg\ 2 dx 
2 Jb Jb 

= I azg(y ' zg ■ v) dx H / |Va| 2 Zgdx / |Vzg| 2 

Jas 2a* J B 2 J B 

< [ azg(Vzg ■ v) dx + /" z ^ d x — ^1 f \\7zg\ 2 dx, 
Job 2a* J B 2 J B 

where ^ is the external normal to dB and in the case n = 1, 3£> = {xi, X2}, £1 < £2 5 
we mean 

az e (Vz e • 1/) dx = a(9, t, x 2 )zg(t, x 2 )zg tX (t, x 2 ) - a(9, t, xi)z e (t, Xi)zg tX (t, xi). 

Substituting (|4~5|) . (|4~7jl in we obtain 

— / —Zg(t,x)dx< — j- f \V zg\ 2 dx + J azgiy 'zg ■ v) dx 
dt J B 2 4 J B J dB 

+ (K B + & + ^)f B 4*. (4.8) 
1 h\\ L ~ in) +C 2 \B\t 2 \\<j\\ L ~ m . 



2 

By Lemma l4.1l and the assumptions on B 



-[ \S7z e \ 2 dx < - 1 / z^dx + -^— I [ z%dx, 
Jb ko\B\ 2/n Jb \B\ 1,n JdB 

so by Lemma EU (Oj) and 



(4.9) 



4Ao|B| z/ "7 7 B 



+ / azgiVzg ■ v) dx H a * , / zi dx (4-10) 

Job A\BV ,n Job 
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< [K 8 + — 7 - + " - " L 5,- / 2| da; 

a, 2a, 4A |B| 2/r V is 

lklU~(TC) + c? I- 8 ! * 2 lklU°°(K) 



+ / a^g ( Vzq ■ v) dx + 



o 



for some constant a' > assuming that, say, e.g., \dB\ < 1. We will eventually 
choose |-B| < £ < 1 sufficiently small (maybe dependent on ||cr||L°°(-R,)) and A 
sufficiently large (independent of ||er|| jw-m) so that 

C\ \B\ e\\a\\ L ^ + < _|_.| k | UooW . (4.11) 

Furthermore, 

K9 2 B Ikll l-(7J) + ^ aze(Vz e • i/) < ^ |S| ||er|| ioc(7J) . 
There exists to > (independent of ||c||l°°(7j)) such that 

"* „ K ? II Vail? oo u) 

Ks-- 1 - n > — 57^ (4-12) 



4A |B| 2/ " 2a, 2 |£| 2/ 

Substituting and gJZ) in jCTty . we have 

_d 



z e (t,x)dx < - / z g {t,x)dx + A\B\ ||<7|| ioo(7i) 



\B\ L/n w \B\ 

By the Gronwall inequality and H3.8[) . we have 



+ iAT + T^TM W - Wli-(R») • ( 4 - 13 ) 



zj(t,x)dx < exp(- 2/n ) / 2 2 (0,a;)cfa: 

+ Aex p(-^pr)/ ex p(^Ar)l B l lklU-(*)^ 

/ wi s f , cjt s At 2 

/ Ldt , f l , LOT \ oi 2 

+ exp ( - 7^) i exp TgiV" — ^l-(-) dT 



< ex P ( - ^at) / M*) - ^^) 2 dT ( 4 - 14 ) 



+ AJ-^— IMU~(K) Il-exp(-^p 
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— \\<P - i>\\ 2 L oo (M . n) \B\ 1/n [ 1 - exp ( - 



w " ^ ' ■ ■ \ " ^ \B\ 2/n ' 



Observe that, 

\B\ 2/n ( , ut ,\ , wf 
— 1 - exp ( ^7-) < t, 1 - exp ( 5- 



and, by (E3 and Remark O 

< iiLo(]|a - 6|U«»( W ) + ||V a ■ / - V x • g||L~(7?.) 
+ ||/« ~ 9u\\l^(tz) + \\ h ~ < f < T, 

for some positive constant Kq, then, from (|4.14|) and since |B| < 1, 

Zg(t,x)dx< / (tp(x) — ij}(xj) 2 dx 



B 

+ K Q A + ij t 2 (||a - fe|| L oc (TC) + ||V S • / - V, • 5 |U=c (TC) 

ck' I -B I ^ n 

+ \\fu -3u|U«>(7i) + l|/i- fc|U«>(w)) H llv - V'IIl-cr") 

+ X A (l + ~J |S| 1/n t 2 (||a - b|| L = (w) + HV, • / - V x • 

+ ||/« ~ 9u\\l<*>(h) + \\ h ~ fc IU°°(7Z))- (4-15) 

Let now E D E be & connected set such that interior of E contains the closure of 
E, dist(dE,dE) > 0, and \E\ — 2 \E\. Since the closure of E is compact, we can 
cover it with finitely many balls B±, . . . , B m C R n , that is, E C Uj-Bj. We may 
choose the balls such that UjBj is contained in the interior of E, and thus 

k 

\U j B j \<Y,\B j \<\E\ = 2\E\. 
j=i 

We assume that both |di3j| < 1 and |L?j| < 5 < 1. Thus the result l|4.15|l holds and 
we may sum the inequality over all balls B\ , . . . , B m C 1" , which yields 



J 4(t,x)dx<c(\E\ + \E\^ n ) \\<p-4,\\ 



2 

L°°(R") 



+ C \E\ 1/n t 2 (\\a - b\\ L ^ (n) + ||Vx • / - Vx • 0|U~ (W) 
+ II/™ - 9u\\l^(k.) + \\h- k\\L^(iz)) (4-16) 
which proves (|4.2() . □ 

This proves the following result. 

Theorem 4.3. Fix T > 0. Lei u = it(t, x) and v = v(t, x) be the classical solution 
of tm and (|2.2|) . respectively, with a = a(t,x,y,q) and b = b(t,x,y,q) satisfying 
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(Hi), f = f(t,x,y), g = g(t,x,y), h = h(t,x,y,q), and k = k(t,x,y,q) satisfy- 
ing (H2), and ip and ip satisfying (Hz). Then there exists a positive constant C 
depending only on T, n, a*, a*, 6», b* , k\, k2, and k 3 such that 

\\u(t, •) - v(t, -)\\ LHE) < C( \E\ 1/2n + \E\ 1/2 )y - ^lU-fRn) 

+ Ct(\\a- b\\ L o° (no ) + ||V X • / - V x ■ gh-pi) (4.17) 

+ Wfu - gu\\L°°(K) + \\h~ fc||z,~(TC )) \E\ 1/2n , 

for allO<t<T with K := [0,T] xEx [-K^Ki], TZ Q := [0,T] x E x [-Ki,Ki] x 
\—K2,Kq\ where E C R n is bounded connected set with Lipschitz boundary. 

Proof. Direct consequence of (|3.2[) , (13. 4|) and Lemmas 14.21 15.11 and 15.21 □ 

5. Estimates in L p (E) 

We want to extend the estimate of Theorem 14. 31 to general p. 

Lemma 5.1 (Case 1 < p < 2). Assume (Hi), (H2), and (Hz)- There exists a 
positive constant C3 depending only on T , n, a*, a*, 6*, b* , k\, ki and k% such that 

\\ze(t, -)\\ LHE) < C 3 ^|(2-p)/ffp)+V2n + | £| i/^ | b _ ^| Uoo(Rn) 

+ C 3 t (||o - 6||l-(Wo) + IIVx • / - V x • (5.1) 

+ ||/« - *IU- W + \\h- k\\ L ^ no) ) , 

for each < t < T , E a W l bounded connected set with Lipschitz boundary, 
< 9 < 1 and 1 < p < 2. 

Proof. By the Holder inequality, 

<|£| 1/<? 'K(V)IIl 9(£ ) (5.2) 



|£| W Uz™(t,x)dx\ 



1/9 



with 



So, by ((521, 



2 

9 : 



P ' 2 -p 

p/2 



||^(t,-)lli.( £ )<l^ (2 " , ' )/2 

<\Ef-^ 2 \\z (t,-)\\l 2[EV 

then, by Lemma T4. 21 

N(t, OIU^) < C 2 (^iM/W+Van + i^iVp) _ ^|| 

+ C 3 t(\\a- b\\ L ^ {K) + || V/ - Vg\\ L <~(ii) 
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+ II/. - fclU- W + \\h - fc|U~(*)) ^l^-rt/^+V^n 5 

This concludes the proof. □ 

Lemma 5.2 (Case p > 2). Assume (7~Ci), (H2), and (W3). There exists a positive 
constant C4 depending only on T, n, a*, a*, 6*, 6*, fci, &2 and k^ such that 

\\ze(t, -)\\ LHE) < C 4 (l + &-^)(\E\ 1/np + \E\ 1/P )\\<p- ^|£V) 

+ C 4 (t + i 2/p )(||a - fclU-TO + ||V» • / - V* ■ gU-m (5.3) 

+ ll/u - ff„|U- w + \\h- k\\ Lo , {n) ) 1/P [E] 1 ^ , 

for each < t < T , E C M" bounded connected set with Lipschitz boundary, 
< 6 < 1 and 2 < p < 00 . 

Proof. Observe that 

\\ze{t,-)\\ P LP{E) = J^z p (t,x)dx 

< INe^Olli^Rn) / zj(t,x)dx 

= lke(*r)llz,~ 2 (R»)ll^(*,-)lll3(j5)- 
Since 2/p, (p — 2)/p < 1, by Lemmas 13.41 and 14.21 we have 

\\ze{t^)\\ LHE) < N(*,0llfcKl|^(«,0||^ 

(c^-^t^ + ib-v-H^) 

+ C 2/p ^ l^l 1 ^ (||o - 6|| i0 o (w) + ||V, ■ / - V x • fflU-^) 

+ ||/u - 5u||i«(7Z.) + - fc|U°=(7?.)) 
4 



< 



1/p- 



< (cf- 2)/p i^ 2 )/p + fc. 

+ C 2 2 /p t 2 ' p \E\ 1/np (||a - + || Vx • / - V x • 0|U. (W) 

\ r/p-i 

+ ll/u - ffulU-(TC) + ||/i - k\\L^(TZ)) 

= (c[ p - 2)/p t(P-V/r + k 4 ^cl /p (\E\ 1/np + \E\ 1/P ) 
x ||^-^|| 2 /f (R „ )+ C 2/p (^- 2)/ ^ + M 2/p ) \E\ 1/np 
x M|o - &|U°°(w) + HVz ■ / - V s • g|U»o(TC) 

+ ||/u -.9u||l-(k) + ||/i - fc|U~(TC) J 



1/p 
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where k& is a positive constant such that 

< ki , 2 < p < oo . 

Since the maps 

2 < p < oo C^- 2)/p , C 2 2/p 
are bounded the proof is done. □ 

The following theorem summarizes the result in Theorem ^31 with the extension 
to general p. 

Theorem 5.3. Fix T > 0. Let u — u(t, x) and v = v(t, x) be the classical solution 
of (I2.1|l and <|2.2|) . respectively, with a = a(t,x,y,q) and b = b(t,x,y,q) satisfying 
(Til), f = f(t,x,y), g = g(t,x,y), h = h(t,x,y,q), and k = k(t,x,y,q) satisfy- 
ing (Ti-2), and if and ip satisfying (W3). Then there exists a positive constant C 
depending only on T, n, a*, a*, &*, b* , k\, ki, and k 3 such that 

\\u(t, •) - v(t, -Jlli-CB) < A E (t)\\<p - i>f L p z im 

+ B(t)(\\a - b\\ L <~ m + HVfl, • / - V x ■ ff|U»( K ) (5.4) 

\E\ r ' p , 

withTZ := [0,T] x E x [-Kt,Ki], TZ := [0,T] x E x [-iTi,^] x [-ifa.ifg]- Here 
h */l<P<2, = f^ + 5 L, z/l< P <2, 

i/2<p<oo, p ' [i, i/2<p<oo, 

( |^|(2- P )/2p+l/2„ + \ E \VP^ tfl< P <2, 

(1 + t lr-*)/p)(\E\ 1/np + \E\ 1/P ), if2<p<oo, 

t, ifl<P<2, 
(t + 0P), if2<p<oo, 

for all < t < T , where E C M. n is bounded connected set with Lipschitz boundary 
and 1 < p < oo. 

Proof. Direct consequence of (I3.2JI , (|3.4(l and Lemmas 14.21 15.11 15.21 □ 
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